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1 BACKGROUND AND OBJECTIVES 


► Two-equation eddy-viscosity models (TEM’s) are the most cost effective 
for the purposes of applied CFD. Give best accuracy vs. cost balance. 

► There is a lot of confusion about true strengths and limitations of TEM's 
especially that of standard k-£ model. 

► We have embarked on extensive study of TEM’s over wide range of flows: 

> Identify true strengths and limitations of standard k~e model. 

> Evaluate other TEM’s. 

> Assess emerging models and novel modeling trends. 

> Identify key areas requiring further research. 

► This talk provides brief review of TEM's from perspective of applied CFD. 

> It provides objective assessment of both well-known and newer models. 

> It compares model predictions from various TEM’s with experiments. 

> It identifies sources of modeling error and gives historical perspective of 
their effects on model performance and assessment. 

> It recommends directions for future research on TEM’s. 


REMARK: 

► Many reported poor predictions of TEMs are primarily due to 
combination of improper choice of near-wall model and over-diffuse 
numerics. 

► TEM performance can be much improved form further research in: 

> Length scale deter minin g equation. 

> Advanced (Anisotxopic/Nonlinear) Eddy-viscosity models. 
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INTRODUCTION AND BACKGROUND 


O About FDI 

► Over 10 years in business. 

► Primary product F3DAP (Fluid Dynamics Analysis Package). 


O About FIDAP 

►* First commercial general-purpose finite element CFD program. 

► Models wide range of flows. 

► Over 700 FIDAP licenses worldwide. 


O FIDAP Turbulence Modeling Capabilities 

► Based on two-equation eddy-viscosity models: 

> Standard &-£ model (Launder and Spalding). 

> Extended k-e model (Chen and Kim). 

> RNG k-z model (Yakhot, Orszag, Thangam, Gatski and Speziale). 

► Low-Re near-wall modeling based on two-layer approach: 

> Viscous sublayers spanned by single layer of specialized elements. 

> van Driest’s model used in viscous sublayers. 

t> Interpolation functions based on universal flow profiles. 

Latest turbulence modeling enhancements (to appear soon): 

> Anisotropic eddy-viscosity models. 
t> Wilcox’s k-<o model. 

> Anisotropic version of the standard k-z model. 


O Typical Industrial User 

► Design engineer. 

► Trained in fluid mechanics and heat/mass transfer. 

Familiar with range of flows of interest to his/her organization. 

► NOT CFD expert. 

► Little or no background in turbulence modeling. 


O Turbulence Modeling Requirements of Applied CFD Codes 

► Optimal balance of cost and accuracy: 

t> Turbulence modeling overhead of critical concern. 

> Overall accuracy of ± 15% adequate for most cases. 

► Consistent performance over wide range of flows: 

> Heat/mass tr ansf er 

> 2-D and 3-D (Cartesian, axisymmetric) 

> Complex geometries 

> Transient flows 

► Adaptable to a wide range of complex flow physics: 
o Low-Re effects 

> Variable density/compressibility effects 

> Combustion 

> Two-phase 

► Minimum level of user input/intervention: 

> No fine tuning model coefficients and/or solution parameters. 

> No physical input other than boundary and/or initial conditions. 

► No geometry dependence: 

> Distance to wall and/or y+ dependence. 

► Stable numerical characteristics. 
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TURBULENCE MODELING CONSIDERATIONS 


O Key Modeling Issues 

1- Accurate modeling of mechanisms governing p u/u' .pu^F, pujc£ . 

a) Pressure-scrambling 

b) Body forces 

c) Transport effects 

d) Dissipation 

2- Accurate modeling of characteristic turbulent length scales. 

3- Accurate modeling of low-Re near-wall phenomena. 

O Optimal Level of Turbulence Model for Applied CFD 

► Second-Moment Closures (DSMC’s) and (ASMC’s) 

(+) DSMC’s ideally suited to modeling aspects l-a,b»c above, however, 
(-) DSMC’s costly, especially in 3-D in presence of heat/mass transfer* 
(-) Geometry dependence in current pressure-scrambling models. 

(-) ASMC’s perform erratically (1-c above not well modeled). 

(-) ASMC’s numerically less stable (stiff equations). 

► Two-Equation Eddy-Viscosity/Diffusivity Models (TEM’s) 

(+) Least costly. 

(+) No geometry dependence (except some low-Re TEM’s). 

(+) Numerically more stable. 

(-) Conventional TEM’s not suitable for modeling effects l-a,l-b,&l-c. 
(+) Room for significant improvement in predicting effects of complex 
strain and anisotropy through the combined use of improved length 
scale equations and advanced eddy- viscosity models. 

(-) Transport effects (1-c), however, cannot be directly predicted. 


LENGTH SCALE DETERMINING EQUATION 


O THE STANDARD fc-e MODEL 
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c ^ = 0.09, Cj = 1.44, c 2 = L92, a k = L0, a £ = 13 

► Remarks on standard k-z model: 

> Use is made of Boussinesq’s “isotropic” viscosity model. 

> Fine scale isotropy is assumed in modeling e equation. 

t> Is high- Re model. Must be used with suitable near-wall sub-model. 

> Many reported poor predictions are due to improper choice of near- wall 
model, mesh density, discretization scheme and boundary conditions. 

t> Model predicts much better than commonly believed, if used properly. 

o It does however have its shortcomings in predicting difficult flows 
involving strong anisotropy and/or non-equilibrium effects - it tends 
to be over-diffuse. It predicts flatter flow profiles, shorter recirculating 
zones, and occasionally does not predict subtle separation bubbles. 
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T KNGTH scale determining equation 

o THE EXTENDED k-e MODEL OF CHEN AND KIM 

► Employs modified e equation containing extra generation term. 

► Rationale is that in addition to turbulence time scale fc/e, there is 
further time scale pkJG characterizing response of £ to mean strain. 
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Cp - 0.09, c, = 1.15, c 2 = 1.9, c 3 = 0.25, o t = 0.75, o t = 1.15 

^ Remarks on extended k- e model of Chen and Kim : 

> Is high-Re turbulence model. Needs near-wall model. 

> Gives similar predictions to standard model in equilibrium flows. 

> We find Chen and Kim's (1987) recommended model produce 
predictions that are too under-diffuse in confined flows. 

> We have tuned constants c, = 1.35 and c, = 0.05 to improve performance 

> Revised model gives better results for some well-known benc hmar k 
flows, but improved predictions over standard model are not realized 
consistently. More experience and possibly fine tuning is needed. 


LENGTH SCALE DETERMINING EQUATION 

O THE RNG k-t MODEL 

► RNG k-z model has undergone two major revisions. 

► Latest version due to Yakhot, Orszag, Thangam, Gatski, and Speziale 



and 


CvL = 0.085, q = L42, c 2 = L68, o t =a e = 0.7179, q 0 = 438, p = 0.015 

Above version is high-Re turbulence model. Needs near-wall model. 

► Most testing of model has been done with simple near-wall model. 

► Our testing of model with more accurate near-wall model indicates 
that RNG model is often under-diffusive in internal flows and can be 
very over-diffusive in some external flows. 

► We have tuned model constants and obtained better overall predictions. 

C>1 = 0.0865, C| = 1.45, c 2 = 1.83, o k = 0.8, a c = 1.15, q 0 = 4.618, p = 0. 17 

► Revised model gives better results for some well-known benchmark 
flows, but improved predictions over standard model are not realized 
consistently. More experience and possibly fine tuning is needed. 
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LENGTH SCALE DETERMINING EQUATION 

O Additional Remarks on RNG k-z Model: 

► Interesting development though no major breakthrough. 

► Most model constants are predicted from RNG theory. 

► In applying RNG theory it is assumed that turbulence field has very 
gide spectaim and that in ertial snb-rep ge is i sotrop ic. 

► Values of model constants predicted by RNG theory are approximate 
owing to simplifying assumptions made in applying RNG method. 

► Model predictions critically dependent on additional term R. 

► The R term reflects proposed contributions from fine sea 1e anisotropy . 

► The R term is not derived and modeled using RNG theory. 

► The R term has essential similarities with extra term in e eq'n 
of extended A-e model of Chen and Kim. 

► Latest model does not predict von Harman constant. 

► The most notable fact about the RNG k-E model of YOTGS is 
that it challenges the notion of fine scale isotropy of turbulence 

> Thus e (and consequently the characteristic turbulent length scale) 
is assumed to be significantly influenced by the fine scale structure. 
These effects are heuristically modeled via the time scale ratio tj. 

> It is interesting to note that the assumption of fine scale anisotropy 
used in modeling R conflicts with notion of a wide and isotropic 
turbulent spectrum used in applying RNG theory to rest of model. 

> It is more likely that the turbulent length scale is influenced strongly 
by large scale anisotropy as characterized by the anisotropy 
tensor a * 

> Anisotropic eddy -viscosity models can provide estimates of a* which 
can be used to design improved length scale determining eqn’s. 


| ADVANCED EDDY-VISCOSITY MODELS (Beyond Boussinesq) 


O Anisotropic Eddy Viscosity Models (AEVM’s) 

► There has been renewed emphasis in developing AEVM’s. 

► Lead to better approximations of the normal and shear stresses and 
therefore turbulence anisotropy effects. 

► In addition to more accurate modeling of , AEVM’s could 
potentially be used to improve modeling of: 

> Length scale determining eq'n. 

t> Generation rate of turbulence energy. 

► Examples of AEVM’s are: 

> Lumely (1970) 

> Spezdale (1987) 

t> Yoshizawa (1984), DIA 

> Rubinstein and Baron (1990), RNG 

> Taulbee (1992) and Speziale (1993), derived from DSMC’s 

> Launder ( 1993) 

► Remarks: 

> Potential of models have been demonstrated using simple tests. 

> Improvements in accuracy often of second-order in magnitude. 

> Not been extensively tested especially for swirling flows. 

> Anisotropic models not yet extended to turbulent heat/mass fluxes. 

> We are presently investigating AEVM’s of Speziale (1987) and 
Launder (1993). 
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THE LOW-RE NEAR-WALL MODEL 


O Wall Function Models 

► Produce over-diffuse solutions in off-equilibrium boundary layers. 

> Often fail to predict separation or vortex shedding. 

► Unfort una tely still in extensive use in applied CFD codes. 

O Specialized Finite Element Model (FIDAP) 

► Is essentially two-layer model. 

► Avoids fine near-wall mesh via use of one layer of specialized elements. 

► Employs van Driest's low-Re mixing-length model in near-wall layer. 

► Combines low cost of wall function models with accuracy of two-layer 
models. 

► y + dependence confined to single layer and transparent to user. 


Remarks: 

► Most of historical testing and verification of TEM*s has been done using 
wall functions. The excess diffusion has lead to much confusion in 
assessing TEM’s. 

► Proper assessment of TEM's requires at least two-layer models. 

► Wall function approach is simply unacceptable for applied CFD. 


IMPACT OF DISCRETIZATION ERROR 


O Sources of Discretization Error: 

► Grid refinement (grid convergence). 

► Location of computational boundaries (e.g., outlet, inlet, entrainment). 

► Choice of discretization scheme in space and time. 


Remarks: 

^ Effect of discretization error hag received less attention in turbulence 
model development and testing. 

► Most serious source of error results from discretization of advection 
terms (i.e., the upwinding scheme). 

► Common but dangerous upwinding strategy is used in many CFD codes: 

> Use accurate unbounded scheme in mean fl ow equations . 

> Use inaccurate numerically diffuse scheme in turbulence equations . 

> Overall scheme is stable but often highly diffusive. 

> Most of development and testing of turbulence models has been done 
using above upwinding strategy. 

t> In our computations we employ the accurate streamline upwind 
(SU) scheme in both mean and turbulence equations. 

> Even more accurate schemes are available which are based on 
Petrov- Galerkin finite element formulations. 

► Accurate schemes must be used in both mean flow and turbulence eo’s. 
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| NUMERICAL RESULTS 

O Free Jets 

► Round jet 

► Plane jet 

O Internal Flows with. Separation 

► Flow past backward facing step (Kim et. al) 

► Flow past step in chann el with diffuser wall (Driver and Seegmiller) 

► Flow in pipe expansion (Szszepura) 

O Transient Flow (V ortex Shedding) 

► Flow past square prism (Lyn) 

O 3-D Flow 

► Flow past passenger car models 
REMARKS: 

► Five sets of model predictions are presented: 

> Standard k-z model 

t> Extended k-z model (original) 

> Extended k-e model (revised) 

> RNG k-z model with (original) 

> RNG k~z model with (revised) 


free jets 

The Submerged Plane and Round Jets 



Plane Jet 

Round Jet 




°fo error 

Experiment 

=0.105 


=0.095 


Standard k-z model 

0.104 

-1 

0.112 

18 

Extended k-e model (original) 

0.10 

“5 

0.10 

5 

Extended k-z model (revised) 

0.102 

-3 

0.104 

9.5 

RNG k-z model (original) 

0.131 

25 

0.157 

65 

RNG k-z model (revised) 


-A 


19 
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TURBULENT FLOW OVER BACKWARD FACING STEP 
Kim et al Test Case: Re = 45000 




% error 

Experiment 

7.0 ±0.5 


Standard k-z model 

6.5 

—7.1 

Extended k-z model (original) 

8.4 

20.0 

Extended k-z model (revised) 

7.1 

1.4 

i ii 1 1 — 

7.5 

7.1 

RNG k-z model (revised) 

7.46 

6.6 



TURBULENT FLOW OVER STEP IN CHANNEL WITH 
DIFFUSER WALL 

Driver and Seeg mille r Test Case: Re = 36000 



Angle 1 


6 degrees 

mm 


mm 

% error 

Experiment 

rrm 


Q3B 


Standard k-z model 

5.3 



-18.5 

Extended k-z model (original) 

6.6 

6.5 

9.55 

17.9 

Extended k-z model (revised) 

5.76 

-7.1 

MEM 

-8.6 

RNG k-z model (original) 


-0.5 

8.33 

2.8 

RNG k-z model (revised) 

6.11 


8.33 

2.8 
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TURBULENT FLOW IN PIPE EXPANSION 
Szczepura Test Case: Re = 890,000 



X* 

% error 

IE 

9.51 


1 Standard k-e model 

9.59 

0.9 


12.44 

30.8 

1 Extended k-e model (revised) | 

10.6 

11.5 


11.35 

19.5 

RNG k-e model (revised) | 

11.39 

20 



TURBULENT FLOW PAST SQUARE PRISM 
Lyn's Test Case: Re = 21400 



Stronhnl No. 

n 

C, 

Experiment 

0.132 ± 0.035 

= 2.0 

NJL 

Standard k-e model 

0.128 

1.68 

0 

l^4.|.U.VJ4..!.!.UIM.U.UU 

0.131 

2.56 

-0.07 

Extended k-e model (revised) 

0.135 

2.014 

0 

RNG k-e model (original) 

0.133 

2.38 

-0.07 

RNG k-e model (revised) 

0.133 

1.9 

0 
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| CONCLUSIONS AND DIRECTIONS FOR FURTHER RESEARCH 


► For applied CFE), TEM's strike balance between accuracy and efficiency. 

► The use of inadequate near-wall models and over-diffiise numerical 
schemes obscures true performance characteristics of TEM's. And this has 
lead to much confusion in evaluation of TEM's. 

► Consequences of using better near- wall model and accurate numerics are: 

> Standard k-e model performs much better than commonly believed. 

> Extended k-e model with original set of model co nstan ts produces 
tinder-diffuse predictions. 

t> RNG k-e model with original set of model constants gives predictions 
that can be both under-diffusive or over-diffusive depending on flow. 

> The extended and RNG models with revised set of model constants 
perform better than with original set of model constants. 

► Newer models are quite promising, but do not yet perform consistently 
better than standard k-e model. 

► Significant advances in TEM capabilities may potentially result from 
further research in two key areas: 

> Advanced constitutive-type laws for the Reynolds stresses: 

• AEVM’s appear to be best candidates. 

> Improved length scale determining equation: 

• Better modeling of off-equilibrium effects. 

• Better modeling of large-scale anisotropy effects. 
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